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Abstract 

A fairly general Lorentz-covariant quark model of mesons is constructed. It has several 
versions whose nonrelativistic limit corresponds to the well-known Isgur, Scora, Grin- 
stein, and Wise model. In the heavy-quark limit, the covariant model naturally and 
automatically produces the heavy-quark symmetry results for meson decay constants 
and semileptonic decay form factors. The meson decay constants and the Isgur- Wise 
functions are calculated for various versions of the covariant model and compared with 
other estimates. A general and adaptable structure of the covariant model ensures that 
it can be used to describe transitions involving light and/or heavy mesons. 
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I Introduction 



The well-established, simple, and often used, nonrelativistic quark model of Isgur, 
Scora, Grinstein, and Wise (ISGW) [1] has also been employed [2,3,4,5] in the inves- 
tigations of heavy-quark symmetry (HQS). Although the ISGW model helped in HQS 
investigations , this nonrelativistic model was not capable [4,5] of properly reproducing 
all of the heavy-quark effective-theory (HQET) relations among semileptonic meson 
decay form factors. It had to be "relativized" to some extent [2,4,5]. Moreover, even 
in the original paper [1] some compensation for relativistic effects had to be introduced 
with meson wave functions. In this way, useful insights in the HQET were gained and 
subleading corrections of order A/tuq were estimated [5] (Here A ~ Aqcd and mg is 
the heavy-quark mass). 

Thus it seemed useful to develop a fully covariant model that, in the nonrelativistic 
limit (NRL), goes into the ISGW model. It turned out that such a covariant model 
can, to a great extent, retain the simplicity, which was an endearing and useful feature 
of the nonrelativistic model [1]. 

The covariant model can have a fairly general form [6] that can be, if wanted, spec- 
ified in such a way as to lead to the ISGW model in the NRL. The given covariant 
formulation allows reasonable freedom in the selection of model parameters and model 
meson wave functions. They can be selected to reproduce a particular Isgur- Wise func- 
tion (IWF), what might provide a good basis for the calculation of A/mg corrections. 

An important feature in all variations of the proposed covariant quark model (CQM) 
is the description of valence quarks (antiquarks). They are parametrized by the on- 
mass-shell Dirac spinors, as it was the case in earlier models [1] and in all subsequent 
usages [2-6]. In a covariant model, such a description might lead to difficulties with 
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the covariant definition of a meson mass M. As shown in the next section, this can 
be resolved by introducing a scalar function that represents the neutral sea with a 
momentum K and vacuum quantum numbers [6]. This is an attractive feature, as sea 
contributions must figure in a description of a hadron. In the present model, which 
takes into account only fluctuations involving valence quarks (antiquarks), the sea is 
described in the simplest possible way, as a physical vacuum. The sea momentum 
function F(K) has a particularly simple form if one wants to define a model that in the 
NRL goes into the ISGW model. 

In the third section of this paper the meson decay constants and the IWF are 
calculated for this version of the CQM inspired by the ISGW model. It turns out that 
the covariant formulation takes care of the relativistic effects, which previously had to 
be compensated for by a phenomenological parameter k (see Fig. 1). The sea function 
in this version of the CQM is just a Dirac delta function, which ensures that a meson 
has a properly defined on-mass-shell four-momentum P, (P 2 = M 2 ). 

However, one could use a nontrivial F(K) function in the CQM. The form of such 
a function would influence the model description of the physical quantities. This is 
illustrated in the concluding sections of this paper by calculating meson decay constants 
and the IWF for a Gaussian F(K). 
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II Relativistic model(s) and the ISGW limit 



A meson H with the four-momentum P and the mass M is covariantly represented 



Here, the index d refers for concreteness, to a light d antiquark, whereas the index Q 
denotes any of heavy quarks. The Dirac functions such as 5{p 2 — nig), combined with 
the corresponding step function 0(e), ensure that valence quarks are on the mass shell. 
This is a characteristic feature of the ISGW model. One can select, if desired, a quark 
wave function which goes into the ISGW wave function in the NRL. Various 

momenta in (2.1) are 



by 




C,S1,S2 



■ d 4 q 5{q 2 - m 2 d ) 0(e) d A K F(K) <5 (4) (p + q + K — P) Q(E) 0(/ ± ) 



(2.1) 



■ u Q, Sl ip)l^d, S2 (q) dj(q, c, s 2 )b^(p, c, si)|0> 



r ± (p) = i 



pv(p ■ i) 

~W 2 



(2.2) 



p» = (e,p) ; g" = (e,g) 



l±(M) = (0, 



P-Q 
2 



The ISGW limit is obtained if the wave function is selected as 
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Wl) = -^e+^WI (2.3) 

Here (5s corresponds to the variational solution [1] with the harmonic-oscillator wave 
functions. 

The sea function F(K) ensures that the meson mass M can be covariantly defined. 
Without F(K), the Dirac delta function 5^(p + q — P) leads in the rest frame (P = 0) 
to the momentum-dependent mass [1] M = (p 2 + ran) 1 / 2 + (p 1 + m 2 d ) 1 ^ 2 . Generally, one 
has the same freedom in selecting F{K) as one had in selecting (f>(l±). However, the 
ISGW state vectors of the weak-binding limit (WBL) [1] will be obtained if a simple 
form is selected: 

pti pv 

F{K) = SP>[K» - ^{^{P ~ (P + ?)),] (2.4) 
In the meson rest frame: 

if°(P = 0) = [M-e-e]=^(p,g) ; K{P = 0) = (2.5) 

Obviously, as /j>k(p,Q) is not always positive, the K does not correspond to a physical, 
on-mass-shell particle. It can be associated with some sea contribution. This contribu- 
tion can, in principle, have a less naive form than (2.4), which has been inspired by the 
ISGW limit. For example, one could try (see Sec. IV) the form 

F(K) = 6W[K» - ^(^(P - (P + ?))„] e~" K2 (2.6) 

It also leads to the ISGW model in the NRL. 

After performing the integrations d 4 K, dp , and dq° in (2.1), one is left with 

\H(E,P,M)> = N £ /^^^ Kp + 9 r-^(^-3£±il)] 
• ^Q, sl (p)^d,sM d + M c, s 2 ) 6j(p, c, Sl )|0 > (2.7) 
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Using the notation 



vl = v"-^- 2 {v-P) ; pf = ^(p.P) (2.8) 

one realizes that the Dirac delta function in (2.7) constrains the orthogonal components 
of the quark four- vectors, i.e., 

Pi + ql = (2.9) 
In the meson rest frame this gives the ISGW relation 

p + q = (2.10) 
By rewriting the complex 5 function in (2.7) one obtains a more manageable form 

\H(E,P,M) >=NJ2[ d^d^ E 1 *< 3 >tf+jf- ^(p||)T) 

■Hl±) ^ Q ,sMl^lsM d+(q,c,s 2 )b+(p,c, Sl )\0 > (2.11) 
Here the quantities T and p\\ are 



1 J m d- m Q +P\\ P-p 
= 1 + 1 _ 1 • P || = nr / 

In the NRL and the WBL [1], 



T=1 + S ' { ^ 



-»2 ^ 2 ->2 ^ ^ 2 

P «m Q ; g « m d 



E ^ M ; e — > m Q ; e — > m d 



Pp P 
P||^ m Q-^ I 9-^ -P+(^Q + ^)^ (2-13) 



Pp 

(Z^) 2 -> 2m Q — -p 2 ; u Q , sl (p)l5V d ,s2(q) -> <$ ai) 



-S2 



one finds 



\H(E, P, M) > NR , WB = N £ / d* P [—^e { - p2+2mQ ^] (2.14) 

— * 

P 

•8 Sl -s a d d(-P+ ( m d + m Q ) — ,c, s 2 )b^(p,c, si) |0 > 

The substitution 



f= * + lf P (2 ' 15) 
leads to the well-known [1-5] ISGW form 

|tf(£, P, M) > W b= JV^/ rfV h^e"^] (2.16) 

C,S U S2 J ^ ' PS 

■K,- S2 4{-p> + ^P, c, 82 )bitf + ^P, C, Sl )|0 > 

The full covariant forms (2.1) or (2.11) lead to fully covariant predictions for meson 
form factors in the CQM, as shown below. These states can be covariantly normalized. 
With 



< H(E',P',M)\H(E,P,M) >= IE <5 (3) [P - P') (2.17) 
one finds for the state (2.7): 



e e 1 \ML(] _ ^12 fElLL^a 



p-q + tym^ p _ P% = 2M 3 ( p _ p (2 . 18) 



[X 1 P(eP-Ep) j m Q m d ' n ?=-p+M(p\\)T 

EM^m 2 d -m 2 Q +p 2 

From this expression, N(P) is 



^) = §^(0) (2-19) 



where N(0) can be calculated numerically. 
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Ill Meson form factors and the Isgur-Wise function 



The form factors for B — > D(D*) semileptonic transitions are defined in the standard 
way: 

< D+(P D )\cYb\B°(P B ) >= j^[U{Pb + PdY + f-(Ps ~ PdY] 



< D*+(e, P D )\cYb\B°(P B ) >= -^g € ^X(PB + Pd) p (Pb ~ P D )« (3.1) 



< D*+(e, P D )\af>y 5 b\B\p B ) >= j^[fe^+a + (e*-P B )(P B +P D r+a^e*-P B )(P B -P D Y] 

Here the vector meson state \D* > is obtained from (2.1) by replacing 75 by (£7). For 
example, one finds that 



< 



D^P D )\cYbmP B )>=3(2nYN D (P D )N B (P B ) £ / d 3 p'd 3 p ^ "'" 

Sl,82,8 1 



e e' 



] ~itj:5 {3) [-p+ ^(Pb\\)T b +P> - ZrMTn] (3.2) 



•0d (Id±)(I>b(Ib±) [-Vd,82 (v'hnUc^ {p')u c ^ {p')^u b)Sl {p)u b<Sl {p)^v d ^ (q')\ 

Here 

E D e' - P D p> E B e - P B p 

PD " = W D PB " = M B 

^m 2 -m 2 c + (p m ) 2 ^Jm 2 d -m 2 + (p m ) 2 
±d — 1 H 1 B — 1 H 
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— * 

In the B-meson rest frame (Pb = 0) this becomes 



< 



Here 



D+(P D )\crb\B°(0) >= 3(2nfN D (P D )N B (0) f d'p'd'p^^^ (3.3) 

J e e e 



E D e>) 



{V»} = ^Tr[ 75 (l + + ) 75 (1 - -£-)] 

8 m c m fe m d 



and 



2 m d mbm c m d m b m c m d m b m c m d m c m b 



*' = eip| 4 (m ' _ (pbii)2)I = ^r CTp[ S ] (3 - 5) 



*° = l^W eXPl 20h <:ml ~ (PD " )2)I 
The trace which determines {V^} is analogous to formula (22) of Ref. [4] and the 4 l m d 
term can be connected with the momentum k of the same reference. That term contains 
the Wigner rotation of the light quark. Comparison with (3.1) allows the extraction 
of the form factors, which can be computed numerically. In order to check Lorentz 
covariance, calculations leading to identical results have also been carried out in the 
D-meson rest frame (Pp = 0). 
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The decays with the vector meson D* + in the final states are described by expressions 
analogous to (3.3). The factor {V^} has to be replaced as follows: 



{V} -> {V^} or {A*^} (3.6) 



Here 



{V - n = T^t^i+^i^ (3.7) 
2m r 2m h 2m r! 



In the heavy-quark limit (HQL), one is tempted to identify the heavy-quark mo- 
menta with the heavy-meson momenta, for example, 



i> = (3.9) 
The Dirac delta-function constraints then determine 



(3.10, 

This means that both valence quarks seem to travel as free particles. Indeed, with (3.9), 
one finds that 



m 2 c -(p m ) 2 = ; <p D = ^ 3/2 (3.11) 

One has failed to account for the Wigner rotation of the light quark [4] and all infor- 
mation on the internal quark momenta is lost. Thus, a more reasonable choice is 



V' -> -tj-Pd + k', \k'\/m c «l (3.12) 

It leads to 



10 



^hql = ^JI^ eX Pl^( m2 c - (PD\\) 2 )} \h Q L 



* 3/ Phql 1 Phql 



J D 

Furthermore, in the HQL, m c — > M D , so that 



k'^p'-P D (3.14) 

An analogous procedure is carried out for the B meson. 

If one had chosen (3.9) instead of (3.12), one would have obtained 

{V} = irr[ 75 (l+ f)-f(l+ f)y B ] (3.15) 
m c m b 

Here the Wigner rotation of the light quark is absent. The expression (3.15) is analogous 
to the expressions employed by Ref. [5]. However, this reference does keep some 
information on the internal quark momenta in the valence-quark wave function, by 
retaining some relativistic terms [5], and thus evades the unacceptable result (3.11). 
Finally, one finds, for example, 

f + (P B = 0) HQL = / d 6 p I[p) 1 + ] (3.16) 

J AM B M D m d (E D + M D )m d 

I(p>) = S(2nfN DHQL (P D )N BHQL (0)(^ 1 ) 2 , <Po hql <Pb hql 

E d E D e'-P D q' 
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Here we have used 



e 



e 



E 



D 



M, 



D 



e' e' I ^ 



g' ^q' ; ► — ; e' = \jm d + q' 

m d m d 



-k' + 



D 



M D \ 



m 2 d +(k> ^— ) 



(3.17) 



m b ' m c M D 



<pB — > 4>B HQL 



-2 

1 <? 



L 2/3 2 



HQL 



Similar expressions are readily obtained for other form factors. Defining 



Fi = f + ; A 1 



1 



M B + M D 

we find the well-known HQS relations [7] 



/; V = (M B + M D )g; A 2 = -(M B + M D )a + (3.18) 



F 1 = V = A 2 



-A, 



(3.19) 



(M B +M D ) 2 

This immediately shows that the definition (3.12) has not introduced any X/itiq cor- 
rections. It only retained internal quark momenta and the Wigner rotation, which is 
necessary if inconsistencies and contradictions are to be avoided [4]. The relations (3.19) 
are valid only in the HQL. Then all form factors contain the same Isgur-Wise function 
(IWF), which is determined by 



£(w) = Rf+ = RF 1 



(3.20) 
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_ 2y/M B M D 
M B + M D 

Here 

, (Pd ■ Pb) 
W — V ■ V = „ ^ „ f — 

M D M B 

— * 

In the B-meson rest frame (Pb = 0) one finds 



Ed \Pd\ 

w 



Vw 2 - 1 (3.21) 



M D ' M D 

The expression (3.20) also satisfies the well-known [7] constraint £(1) = 1 . In Fig.l our 
IWF is compared with that calculated by Amundson [5], who obtained 



m d 



t(w) A = exp[- * (w - 1)] (3.22) 
ZK Phql 

Our curve (solid line in Fig.l) is calculated using the parameters of Ref. [5] 



j3 HQL = 0.42 GeV ; m d = 0.33 GeV (3.23) 

A meson decay constant fn is determined by the expression 
1 



(27T) 



P"fH =< 0| : ^(0)7 m 75^q(0) : \H(E,P,M) > 



= 3N(P) J d P ^ Ee _ p J (h)[ mdmQ W^p^t (3-24) 
Here ^q(x) are valence quark fields [1]. The decay constant f H is easily calculated in 

— * 

the frame P = 0. In order to check covariance numerically, it has also been calculated 
for several different P values. The results have always been identical. 
With the parameters [1,5] 



(3 S(U - C) = 0.39 GeV m u = md = 0.33 GeV 



1.645 GeV (3.25) 



Ps - = 0.42 GeV m b = 4.983 GeV 

(ub) 
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one finds that 



f D = 258, 8 MeV ; f B = 151, 9 MeV (3.26) 
In the HQL, the expression (3.24) takes the form 

1 r>it £ q at <B\[*< M \2 Mm d , n M m d P» + Mq 



In the HQL, one uses the average meson masses 



f^HQL = 3N HQL (P) J d p(_) — -^W^t ] (3.27) 



]g J = 3Mj -+ Mg = 1.975 GeV ; M B = 3M ** + ^ = 5.313 G eV (3.28) 



and (3.23) in order to find the HQS result 



h 

The numerical values 



HQL 



/ D „ 0I = 235, 8 MeV ; / B „ 0I = 143, 8 MeV (3.30) 

which have been obtained using (3.27), are quite close to the result (3.26), showing that 
the model-determined corrections to the HQL are about 5 -r 6%. 
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IV Gaussian sea 



It might be useful to demonstrate the flexibility of the expression (2.1) by selecting 
a sea function F(K) that would be different from the naive choice (2.4). In principle, 
this could be based on some QCD modelling of the sea contribution. However, for 
illustrative purposes a simple example can be selected, which in the NRL and WBL 
goes into the ISGW state. Yet, it leads to noticeably different results when used in the 
CQM. This is 

F(K) = 8V>[K» - ^\^{P - (P + q)) v )\ ■ e- aK2 (4.1) 
A simple arbitrary choice for the parameter a is 

a = 2- 1 (3~ 2 (4.2) 
The integration over the sea momentum K gives 

F = ^|M_*liLM] (4 .3) 

In the meson rest frame (P M = (M, 0)), this goes into (2.5). In the WBL, one finds 
that 



K = ; K° = M — [m d + mg) (4.4) 

Thus the K 2 , dependence has disappeared and (2.1) in the NRL- WBL is again the 
ISGW state (2.16). This conclusion is valid for any meson frame, i.e., any P M . 

In the CQM, all manipulations are exactly analogous to those presented in the 
second section of this paper. In all formulas one has to make the substitution 
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<j>(l±) ^<Pk = <P(h)e- aK2 (4.5) 
In the HQL, one can use (4.3) and (3.12)-(3.17) in order to obtain 



K 2 = {M _ pil _ qil y^q2 

(Mhql -> 0(/± We" ag 1 (4.6) 

Again, all formal deductions (3.17) are repeated with insertions of the factors exp[— aqj] 
associated with the product <t>D HQL • (pB HQL - In the frame Pb = 0, for example, one has 

<Pd hql <Pb hql -> <f>D HQL </>B HQL ■ exp{-a[(q Dll ) 2 + e 2 ]} (4.7) 
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V Numerical examples and discussion 



It is well known that various relativistic wave functions (states) can lead to the same 
state in the NRL. This has been illustrated in Sees. II and IV for two slightly different 
versions of the CQM. However, the different versions of the CQM lead to somewhat 
different estimates of physical quantities. These differences persist even in the HQL. 

Owing to the relations (2.4), (4.1), or similar, the interplay of the sea and the 
valence-quark contribution leads to a meson mass M which is not just a sum of the 
valence-quark masses. However, the values of the quark masses are interconnected with 
the quark wave function (2.3). For any mass change, the variational procedure which 
leads to (2.3) has to be repeated. Or alternatively, the meson wave function <fi can be 
determined in some other model, for example, in a model inspired by the Bethe-Salpeter 
equation. For illustrative purposes, the mass has been changed in the model versions 
determined by (2.3) and (2.4) or (4.1). Such a procedure, admittedly inconsistent, has 
been used just to illustrate the flexibility of the CQM's. 

Table I 

In Table I the meson decay form factors are calculated using (2.4) and the param- 
eters (3.23), (3.25), and (3.28). The spectator antiquark mass m d has been arbitrarily 
changed, as discussed above. As expected, the fn values change with ma, but not 
dramatically, mostly by less than 20%. Deviations from the HQL relation (3.29) are 
more interesting. They amount to abouth 9% when the full CQM is employed. This 
indicates that, in this version of the model, the HQS relation (3.29) presents a very 
good approximation. Absolute values Jb{CQM) are closer (less than 5% difference) to 
f^HQLi as it should be expected with rrib/m c = 3.3. 

Table II 
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The results presented in Table II show that the sea contributions can be very im- 
portant. Whereas the general pattern is similar to that displayed in Table I, the ab- 
solute values of the meson decay constants ju are much smaller. For example, with 
m d = 0.33GeV, / D (4.1)// D (2.4) = 0.59, and / B (4.1)// B (2.4) = 0.70. However, this has 
to be taken more as an illustration of the model flexibility than as a serious prediction. 
The sea descriptions (2.4) or (4.1) are very crude and one should better refer to them 
as to "mock sea" functions. 

Other calculations of the meson decay form factors lead to a broad range of values. A 
relativistic quark model with centrally confined quarks [8] gave values that were smaller 
up to 50% (f D = 130.6 MeV, f B = 90.9 MeV) than our (3.26) values or the values in 
Table I. These values are much closer to our values shown in Table II, which illustrates 
how model-based predictions depend on details of the model construction. The larger f D 
and f B values, as obtained in the present CQM, are closer to the results based on QCD 
sum rules, lattice calculations, and semilocal parton-hadron duality [9] . The estimates 
in the QCD sum rules [10] gave f D = 120 -=- 250 MeV and / B = 90 4 200 MeV. The 
predictions of the lattice calculations [11] are in a similar range fo = 170 4- 230 MeV 
and f B = 140 4 220 MeV. Lower values, f D 80 MeV, f B 130 MeV, remarkably 
close to the values of Ref. [8], were found in a potential quark model [12]. 

The IWF calculated in the CQM defined by (2.4) (see Fig. 1) shows similar behavior 
to the ISGW-Amudson [5] result obtained with the correction factor k = 0.6. Its slope 
p defined by 

£M = l-p>-l)+0((^-l) 2 ) (5.1) 

is about 25% larger than the result of Refs. [5]. This can be attributed to the fully 
relativistic character of the CQM, including Wigner rotation. As shown in Ref. [4] the 
Wigner rotation increases the slope of IWF by about 20%. 

Figure 2 shows the influence of the "mock sea" contribution. The solid curve and 
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the dotted curve correspond to the sea functions (2.4) and (4.1), respectively. The use 
of (4.1) increases the slope p by 10%. One is tempted to assume that the relativistic 
effects (see Fig.l) might play a larger role in the calculation of the IWF than the sea 
effects. However, the model is too crude for such far-reaching conclusions. 

An arbitrary d-antiquark mass change, using = 0.1 GeV instead of (3.23), pro- 
duces virtually the same £(w) curve with either the (2.4) or the (4.1) sea description. 
The slope p = 1.07 is 9% smaller than the result based on (3.23) and (2.4) (see also 
Fig. 1). 

A plethora of p values can be found in the literature. By fitting the data on B — > Dlv 
and using different ansatze for the IWF [4,13-15], the p values have been found to be 
in the range p = 0.92 4 1.57. This is slightly larger than the value based [2,3,16-18] on 
ISGW model [1], p ~ 0.8. Relativistic quark models [8,12] gave p = 1.25 and p = 1.1, 
respectively. The QCD sum-rule estimates [19-26] p = 1.0 4 1.14, as well as the lattice 
computations [27-29] p = 0.714-1.35, are more or less in the same range as the estimates 
[4,13-15] based on the B — > Dlv data. The values of the slope parameter p shown in 
Figs. 1 and 2, which roughly span the range p = 0.93 4 1.29, are similar to fitting-data 
estimates [4,13-15], or QCD sum-rule estimates [19-26], or lattice [27-29] estimates. The 
value obtained in the simplest CQM version of the ISGW model is p = 1.17. 

Our results obtained using (2.4) are connected with the Close and Wambach [4] 
deductions. In their case, the ISGW model was relativized sufficiently to make it 
covariant in the HQL. They have also taken care to include the Wigner rotation of light 
quarks. With ttiq 1 corrections, their approach might also approximately, and adequatly, 
describe lighter mesons (K). 

When one starts with a covariant description, the HQL follows automatically by 
the m~Q expansion. The HQS results are readily obtained in this limit. The Wigner 
rotation is also automatically included in a covariant procedure. Furthermore, the CQM 
can be used for the description of the light mesons or the heavy-light meson transitions. 
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A very general structure of the model, including the sea function F(K), provides for 
its great adaptability and ability to model various physical situations. 
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Table I. Meson decay constants corresponding to formula (2.4). 



m d (GeV) 


f D (MeV) 


f DHQL (MeV) 


f B (MeV) 


f BHQL (MeV) 


0.33 


258.8 


235.8 


151.9 


143.8 


0.3 


251.4 


233.2 


149.3 


142.2 


0.2 


227.1 


223.6 


140.3 


136.3 


0.1 


204.3 


212.6 


131.2 


129.6 


0.01 


187.1 


202.2 


123.6 


123.3 
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Table II. Meson decay constants corresponding to formula (4.1). 



m d (GeV) 


f D (MeV) 


fn HQL {MeV) 


f B (MeV) 


f BHQL (MeV) 


0.33 


157.4 


146.2 


105.7 


89.1 


0.3 


153.8 


144.5 


104.4 


88.1 


0.2 


141.0 


137.8 


99.6 


84.0 


0.1 


127.2 


129.2 


93.5 


78.8 


0.01 


114.5 


120.5 


87.1 


73.5 
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Figure captions 



Fig.l The IWF's are shown as functions of w. The dashed curve shows the IWF 
calculated in Ref. [5] with k = 1. The dash-dotted curve shows the result of Ref. 
[5] with k = 0.6. The corresponding slope parameter is p ~ 0.93. The solid curve is 
obtained using our formula (3.20). Its slope parameter p pa 1.17 corresponds to the 
ansatz (5.1). All IWF's were calculated using the parameters (3.23). 

Fig. 2 Several IWF's are shown. The dashed curve, with p pa 1.07, is obtained by 
using either (3.20) or (3.20) plus (4.7) substitution with m d = 0.1 GeV. The solid 
curve, p ^ 1.17, is obtained using (3.20) and m d = 0.33 GeV, the same as for the solid 
curve in Fig.l. The dotted curve is obtained with the Gaussian sea (4.7). Its slope 
parameter is p pa 1.29. All slope parameters corresponds to the ansatz (5.1). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9509332vl 



